Periodic orbits and semiclassical form 
factor in barrier billiards 



O. Giraud 

Laboratoire de Physique theorique, UMR 5152 du CNRS 
Universite Paul Sabatier 
31062 TOULOUSE CEDEX 04 
France 

February 6, 2008 

Abstract 

Using heuristic arguments based on the trace formulas, we analyti- 
cally calculate the semiclassical two-point correlation form factor for a 
family of rectangular billiards with a barrier of height irrational with 
respect to the side of the billiard and located at any rational position 
p/q from the side. To do this, we first obtain the asymptotic density 
of lengths for each family of periodic orbits by a Siegel-Veech formula. 
The result ^(0) = 1/2 + 1/g obtained for these pseudo-integrable, 
non-Veech billiards is different but not far from the value of 1/2 ex- 
pected for semi-Poisson statistics and from values of -ft^O) obtained 
previously in the case of Veech billiards. 



1 



1 Introduction 



Quantum billiards, that is, closed compact domains in the two-dimensional 
Euclidean plane, are the simplest model of a quantum system corresponding 
to physical instances such as quantum dots or microstructures. The sta- 
tistical properties of the quantum energy levels of such systems have been 
investigated, and it turns out that the statistical quantum behaviour can be 
related to the classical properties of the system. It is believed that systems 
whose classical motion is chaotic have energy levels behaving like eigenvalues 
of random matrix ensembles [8] , whereas the energy levels of systems whose 
classical motion is integrable are Poisson distributed, i.e. they behave like 
independent uniformly distributed random variables [7] . Both numerical ev- 
idence and some analytical results support these conjectures [2, 1, 24]. 
Among systems which are classically neither chaotic nor integrable, some 
systems have been found to display an eigenvalue statistics which is inter- 
mediate between the Poisson and the Random matrix distribution. The 
characteristics of such intermediate statistics are [10] level repulsion, expo- 
nential decrease of the nearest-neighbour spacing distribution at infinity and 
linear asymptotic behaviour of the number variance (which is related to a 
non- vanishing form factor at small arguments). The form factor at the ori- 
gin is equal to 1 for classically integrable systems, to for chaotic systems, 
and it is found numerically to take values between and 1 for intermediate 
statistics, the case ^(0) = 1/2 corresponding to semi-Poisson statistics [10]. 
Numerous quantum systems have been found to display numerically inter- 
mediate statistics: for example, pseudo-integrable systems such as rational 
polygonal billiards (polygons in which all angles are commensurate with 7r) 
[13], or quantum maps [17]. 

An analytical approach to the study of level statistics is the semiclassical 
trace formula, which gives an expansion of the density of energy levels as 
a sum over periodic orbits [3, 20], or families of periodic orbits in the case 
of integrable systems [6]. For diffractive systems, the trace formula can be 
modified to include diffractive orbits contributions [26, 11]. It can be argued 
however (see [12] for a discussion) that only the periodic orbits contribute 
to the semiclassical form factor at small arguments, K 2 (0). The calculation 
of this quantity therefore only requires to find the periodic orbits and the 
areas occupied by the pencils of periodic orbits in a given system. Unfor- 
tunately, this is not a simple task. For instance is is not known whether 
any acute triangle has a periodic orbit. In the case of rational polygonal 
billiards, it has been shown [25] that the number N{L) of periodic orbits of 
length less than L is quadratically bounded, namely there exist C\ and c 2 
such that ciL 2 < Af(L) < C2L 2 , but even for general rational polygonal bil- 
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liards exact asymptotics is not known. There exist however certain specific 
rational polygonal billiards for which more precise statements are known. 
For instance for Veech billiards [27, 28], a special class of rational polygonal 
billiards (whose stabilizer is a discrete cofinite subgroup of SX(2,R)), pre- 
cise asymptotics for M{L) is known, and in [12] it was possible to calculate 
analytically the form factor at the origin for triangular Veech billiards. 

This paper presents the calculation of the semiclassical form factor at 
the origin for a billiard which does not have this special Veech property, the 
barrier billiard. The barrier billiard is one of the simplest pseudo-integrable 
billiards. It was introduced by Hannay and McCraw [21] and consists of a 
rectangle [0, a] x [0, b] containing a barrier described by the segment {e a} x 
[0, ab] with < e , a < 1 (see Figure 1 left). It is a rational polygonal billiard 
with six angles equal to n/2 and one angle equal to 2n. It is therefore a 
pseudo-integrable billiard [5], and the movement in phase space takes place 
on a surface of genus 2. When the height of the barrier is such that a e Q 
then the barrier billiard is a Veech billiard. But when a is irrational the 
billiard loses this property. Nevertheless, from results obtained in [15], it 
is still possible to work out the distribution of the periodic orbits in this 
latter case, and thus calculate analytically the semiclassical form factor at the 
origin, provided the position of the barrier is a rational number with respect 
to the size of the side: e = p/q with p, q G N coprime. We will first devise 
a method to obtain a complete characterization of the periodic orbit pencils 
in the non- Veech barrier billiard (Section 2). We then rigourously derive 
asymptotics for each family of periodic orbit pencils (Section 3), then use this 
result to calculate the semiclassical form factor at small arguments (Section 
4). Previously obtained analytical results show that the semi-classical form 
factor at the origin takes non- universal values between and I. For Veech 
triangular billiards with angles (tt/2, ir/n, n/2 — ir/n), the value /^(O) = 
|(n + e{n))/{n — 2) with e(n) = 0, 2 or 6 was found [12]. For a rectangular 
billiard perturbed by an Aharonov-Bohm flux line, we obtained ^(0) = 
1 — kq. + 4a 2 where a G [0, 1/2 [ is the strength of the magnetic flux and 
k a rational depending on the position of the flux line in the billiard (for 
irrational positions, k — 3) [12]. For a circular billiard perturbed by an 
Aharonov-Bohm flux line, a similar result -fG(O) = 1 — K~a(l — a), with 
k G [0,2] an explicit function of the position of the flux, was derived [18]. 
In the case of the barrier billiard, we obtain 1^2 (0) = 1/2 + 1/q. This value 
depends on the position of the barrier inside the rectangle, which reflects the 
fact that the structure and the properties of periodic orbits strongly depend 
on it. This analytical expression for 1^(0) extends previous results to the 
case of non- Veech polygonal billiards. 
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2 Periodic orbits in the barrier billiard 



The aim of this section is to characterize periodic orbits in a barrier billiard. 
We first begin by the simple case of a rectangular billiard. 

2.1 Periodic orbits in the rectangular billiard 

Let us consider a rectangle of area A = a x b with Dirichlet boundary con- 
ditions. It is easy to work out the density of the lengths of periodic orbits. 
Any orbit in the rectangle can be unfolded into a straight line in a torus (a 
rectangle with periodic boundary conditions) of size 2a x 26; a periodic orbit 
is therefore defined by two integers M and N and has length 

l p = v / (2Ma) 2 + (2Nb) 2 . (1) 

If we restrict ourselves to (M, N) in the upper right quadrant, each family 
of periodic orbits occupies an area AA (2A for the orbit itself, 2A for its 
time-reverse). The number Af(l) of pencils of length less than I is just the 
number of lattice points (2Mo, 2Nb) within a (quarter of a) disk of radius 
I. It has the asymptotic expression Af(l) ~ nl 2 /lQA. The corresponding 
density of periodic orbits is the derivative of J\f(l): 

Pd) ~ £ (2) 
The density of primitive periodic orbits is given by (see e.g. [12]) 

PM ~ £4- (3) 

We want to obtain a similar result for the barrier billiard. In the rest of this 
section we investigate the periodic orbits of the barrier billiard, and Section 
3 leads to Equation (26) which gives the density of primitive periodic orbits 
for the barrier billiard. 



2.2 The translation surface 

Instead of studying directly the barrier billiard itself, we will consider the 
equivalent problem of studying the translation surface associated to this bil- 
liard [19]. A construction due to Zemlyakov and Katok [30] shows that 
the translation surface associated to a generic rational polygonal billiard is 
obtained by unfolding the polygon with respect to each of its sides, which 
means gluing to the initial polygon its images by reflexion with respect to 
each of its sides and repeating the operation. If the angles of the polygon are 
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Figure 1: The barrier billiard and its translation surface 

cti = Tirrii/ni and iV is the least common multiple of the rii, then 2N copies 
of the initial billiard are needed. Here all the angles are multiples of 7r/2, 
therefore only 4 copies are needed, and the translation surface S obtained 
by this construction is represented in Figure 1 (right). In this surface, all 
opposite sides are identified. Any trajectory in the barrier billiard can be 
unfolded to a straight line on the translation surface. The surface S is of 
genus 2: there are two singular angles of measure 4n that we will represent 
respectively by z\ (a dot in Figure 1) and z 2 (a cross in Figure 1). The two 
singularities are traditionally called saddles [16] and a geodesic joining them 
is called a saddle-connexion. 

2.3 Periodic orbits in the barrier billiard 

In this subsection, our aim is to describe qualitatively the periodic orbits in 
the barrier billiard in a given direction. On translation surfaces the peri- 
odic orbits occur in pencils, or cylinders, of periodic orbits of same length. 
These cylinders are bounded by saddle-connexions and are characterized by 
their length and their height. Let us consider a 'rational direction' on the 
translation surface S: 

v = (2Ma/q,2Nb), (4) 
with M and N two coprime positive integers. The length of the vector v is 

l p = v/(2aM/g) 2 + (2WV) 2 . (5) 

Let us label by the integers k — 0, 1,..., q — 1 the positions on the translation 
surface such that the barrier on the " left" of the translation surface in Figure 
1 be at position p — 1 and the barrier on the "right" in Figure 1 be at position 
q — p (see Figure 2). Since the opposite sides on the translation surface 
are identified, then when a trajectory hits the barrier at position p — 1 it 
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reappears at position q — p, and vice-versa. The translation by vector v 
induces a permutation a v of the positions {0, 1, q — 1}. Let us define 

Wi = min {k E N; cr*(P _ 1) ^ {p - 1, g - p}} (6) 

and in the same way 

w 2 = min {k E ~N; a^iq — p) E {p — I, q — p}} . (7) 

A translation by the vector w±v takes z\ to itself and defines a saddle- 
connexion of length W\l p . The second saddle-connexion joining z\ to itself 
starts at position q — p and its length is w 2 l v . 

Figure 2 shows, as an example, the two saddle-connexions going from z\ to 
itself in the direction (9,2) for p/q = 1/3. The translation by the vector v 
induces the permutation (012) i— > (021). One of the saddle-connexions goes 
from the position to itself and has a length l p ; the other goes from position 
2 to itself via position 1 and has a length 2l p . In any direction, there are al- 
ways two saddle-connexions going from z\ to itself, and, in the same way, two 
from z 2 to itself. These four saddle-connexions form the boundary of three 
cylinders of periodic orbits (see Figure 3). The lengths of these cylinders are 




Figure 2: Starting from points of abscissa 0, 1 or 2 ( for q = 3) in the direction 
(M = 9,N = 2), one arrives at 0, 2 or 1: there are two saddle-connexions 
-> and 2 -> 1 -> 2. 

necessarily of the form w±l p , u^p and (w± + W2)l p , with u>j G N, and their 
heights (2b/M)hi are such that h ± + h 3 E Z, h 2 + h 3 E Z and /i 3 — (j5 2 G Z 
for some a = ±1 and 5 2 = {Ma}, the fractional part of Ma. For instance 
in Figure 3, there is one cylinder immediately above the saddle-connexion 
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0^0, one immediately above the saddle-connexion 2 — > 1 — > 2, and the 
third cylinder is below both. 




Figure 3: Three cylinders of periodic orbits bounded by four saddle- 
connexions in the case M = 4, N = 1. 

The results of this section can be summed up as follows. We set 

s 1 = h 1 + h 3 

s 2 = h 2 + h 3 , (8) 

so that Si and s 2 are integers. Then in each direction v defined by (M, N) 
with M and N coprime, there are three cylinders of periodic orbits of lengths 
Wilp and heights (2b/M)hi with i = 1,2,3. The cylinders can be described 
by the following five characteristic numbers: 

- the integers w\ and w 2 (giving the lengths of the two short cylinders 
and the length (wi + w 2 )l p of the long cylinder) 

- the real number h 3 (giving the height (2b/M)h 3 of the long cylinder) 

- the integers s± and s 2 (giving the heights {2b/M)(s\ — h 3 ) and 
(2b/M)(s 2 — h 3 ) of the short cylinders). 

Note that by definition of the s« we need to have < h 3 < min(si, s 2 ). Also 
note that the condition that the sum of the areas of the cylinders be 4A can 
be expressed as SiWi + s 2 w 2 = q. 
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3 Asymptotics for the periodic orbit lengths 



Let us define T as the set of all 4-uples (wi, W2, si, s 2 ) € (N*) such that 
(si, S2) are coprime and S1W1 + S2W2 = g. We say that a direction v belongs 
to the family / e J 7 if the three cylinders in the direction v have the charac- 
teristic numbers wi,W2, si, $2- The goal of this section is to calculate, for a 
fixed family f <E J 7 and a fixed interval I C [0, min(si, S2)[, the asymptotics 
for the number J\fffi(l) of directions v belonging to the family /, such that 
||v|| < I and such that the height h 3 of the third cylinder in the direction v 
belongs to the interval /. 



3.1 Counting periodic orbits 

The asymptotics for the number J\f( g \l) of cylinders of length less than I 
have been calculated in [15]. These asymptotics are obtained by applying a 
Siegel-Veech formula to the space A4 q (l,l) of g-fold coverings of the torus 
with two branch points and area I. If V(S) is the set of vectors associated 
with cylinders of periodic orbits on a 'stable' g-fold torus cover S, then it 
is shown that there is a constant k(S) depending only on the connected 
component A4(S) of M. q {\, 1) containing S, such that 

\V(S) fl B{T)\ ~ ttkT 2 , (9) 

where |.| denotes the cardinal of a set and B(T) is the ball of radius T 
centered at the origin. The constant k is given by the following Siegel-Veech 
formula: for any continuous compactly supported ip : R 2 — > R, 

if dfx = k 1 ip, (10) 



'M(S) 

where (p is the Siegel-Veech transform of ip defined by 

m = Yl vt?) ( n ) 

vev(S) 

and fx is the measure on Ai q {l,l) (Theorem 2.4 of [15]). This Theorem 
applies to the translation surface constructed from the barrier billiard, pro- 
vided S be a stable 5-fold torus cover, which is true only if the height a 
of the barrier is irrational. It is shown that in this case M(S) is the set 
V q (l, 1) C M q (l,l) of primitive torus covers. The following asymptotics are 
then obtained (here we have a factor 1/16 differing from the factor in [15] 
because of our conventions for the counting of the time-reverse partner of a 
periodic orbit): 

^»(0~c^. (12) 
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The constant c is given by 
c=-^-V/i(r) V u 1 u 2 (u 1 + u 2 )mm(s 1 ,s 2 ) ( ^ + ^ + 7 — 

r\q (s 1 ,a 2 )=i x v 7 

si«l+S2"2=l/'' 

(13) 

(the gcd of s, s' will be noted either gcd(s, s') or simply (s, s')), and 

N q = ^2K r ) r ' 2 E u i u 2(ui + u 2 )mm(s 1 ,s 2 ) (14) 

r\q (si>«2) =1 

s 1 u 1 +s 2 u 2 =q/r 

(Proposition 4.14 of [15]). The constant N q is the number of primitive covers 
of degree q of a surface of genus 2 with 2 branch points. 



3.2 Siegel-Veech formula 

The proof leading to Equation (12) can be adapted to any subset of V(S) pro- 
vided it varies linearly under SL(2, R) action, i.e. provided the subset verifies 
V# G SX(2,R), V(gS) = gV(S) (see Section 2 of [14] for more detail). To ob- 
tain the asymptotics for a fixed pair F = (/, I) with / = (wi, w 2 , s±, s 2 ) G T 
and / an interval, / C [0,min(si, s 2 )[, let us define Vf(S) the set of vectors 
v G R 2 defined by (4), such that the triple of cylinders in the direction v 
belongs to the family /, with h 3 G I. Then along the same lines of the proof 
of Theorem 2.4 in [15], one can show that when the height a of the barrier 
is irrational, the translation surface S of the barrier billiard is a stable g-fold 
torus cover and 

|Vp(S') l~l B(T)\ ~ irKpT 2 , (15) 
where the constant Kp is given by the Siegel-Veech formula 

hiv \ Tvs [ <Pf dfi = n F <p F , (16) 
M^U) L )) Jv g (i,i) Jr 2 

with ifip the Siegel-Veech transform 

vev F (s) 

for some continuous compactly supported ip : R 2 — > R. 



3.3 Asymptotics for a family of periodic orbits 

Following the steps leading from the Siegel-Veech formula (9)-(ll) to the 
asymptotics (12)-(13) in [15], we can now derive asymptotics for the num- 
ber of cylinders in each family (/, I) . Recall that Af^j (I) is the number 
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of directions v belonging to a family characterized by the numbers / = 
(wi, w 2 , si, s 2 ), with a height h 3 e /, and such that l p given by Equation 
(5) is less than /. (Note that ftffj(l) is a number of directions and not a 
number of cylinders.) Let Pf,i(l) be the corresponding density. According to 
Equation (15), J\fj q }(l) is proportional to I 2 ; we define the constant Cfj by 

The proof leading to the asymptotics for Afj q )(l) is essentially the same 
as the proof in [15], section 4.4, provided we replace the counting functions 
of the cylinders in [15] by counting functions of directions in which the cylin- 
ders belong to the family (/, I) we are interested in. We take to be the 
characteristic function of a disc of radius e in R 2 . Therefore its Siegel-Veech 
transform tp, as defined by (17), counts the number of directions on S in 
which the cylinders belong to family (/, I) and such that l p < e. For e small 
enough, the Siegel-Veech formula (16) is equivalent to 

= C(2) * f tp F djl, (19) 

where ( is the Riemann Zeta function and 

{1 if the cylinders in the horizontal direction belong 
to the family /, if hs e / and if ||v|| < e 
otherwise. 

(20) 

We define Xf,i ■ ^ 2 i-> R by 

{1 if the cylinders in the horizontal direction belong 
to the family /, if h 3 G / and if ||v|| < e^fq (21) 
otherwise. 

Following [15], we parametrize V q (l, 1) and perform the integration in (19). 
Part of it can be related to the integral over Xf,i, which is f R2 Xf,i( v ) dv = 
ne 2 q. The integration yields 

/ <p F djl=^- V ^Iw 1 w 2 (w 1 + w 2 )\I\. (22) 

qy ' ' r\(wi,W2) 

From [15] we get p,(V q (l, 1)) = N q /q, with N q given by (14). Equation (19) 
finally gives 

c f' I = W —w 1 w 2 (w 1 + w 2 )\I\, (23) 

q r\(w 1 ,w 2 ) 
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where |/| is the length of the interval /. 

Equation (23) shows that Cfj depends on I only through its length. Is 
is therefore convenient to introduce the density of directions p = (M, N) 
corresponding to a family / and such that l p < I and h 3 = h: 

^S(0~c />fc ^, (24) 

with Cf t h given by 

c f* = wT, ^iM^i + W2)e,(h) (25) 

q r\ q 

for any family / = (wi,W2, Si, S2) of T and ftei. The function Of is the 
characteristic function of the interval [0, min(si, s 2 ) [■ The density of primitive 
periodic orbit lengths for the family / G T and h E R is 

3/ 

(26) 

It is easy to verify that the expression (24) of Mjl (/) is consistent with the 
total number M^ q \l) of pencils of periodic orbits with length less than /. This 
comes from the fact that any pencil of periodic orbits contributing to M^ q \l) 
belongs to a certain family / and has a length Wil p < I, which implies that 
l p < l/wi. Therefore 

AT® (I) = W dhj^N${l/ Wl ) 
fer J i=i 

~ e/* c «.si(4 + 4 + kt^)' (27) 

Using Equation (25), we obtain, after integration over h, 

q («1.S2) =1 r\{w\,W2) 
s 1 w 1 +s 2 w 2 =q 



xmm( Sl ,s 2 ) ( — + — + - ■ - . (28) 

\wf w$ {w 1 + w 2 yj 

Making the substitution Wi = rui and inverting the two sums, we get exactly 
the expression given by Equations (12) and (13). 
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4 Calculation of the form factor at r = 



4.1 Definitions 

The spectrum {E n ,n 6 N} of a quantum billiard can be described by the 
density 

d{E) = Y,KE-E n ). (29) 

n 

The two-point correlation form factor is defined as the Fourier transform of 
the two-point correlation function of the density of states: 

/°° He 
-j{d(E + e/2)d(E - e/2)) c e 2 ^. (30) 

Here the product of the densities is averaged over an energy window of width 
AE > 1/d centered around E = k 2 and such that AE «£. If A is the area 
of the billiard, d = A/47T is the non-oscillating part of the density of states. 
The subscript c means that one only considers the connected part of the 
correlation function. It can be argued that in the case of pseudo-integrable 
systems, the leading term of the semiclassical expansion of K 2 (t) at small 
argument (r — > 0) is given in the diagonal approximation by the contribution 
of periodic orbits only: K 2 (t) = K dia,g (r) + O(r), with 

v y 

(see [12] for the derivation of this expression, based on heuristic arguments). 
The sum is performed over all pencils of periodic orbits p of length l p . In 
general, there can be several pencils having exactly the same length: in 
Equation (31), S p is the sum of the areas occupied by all pencils having, when 
(possibly) multiply repeated, a length l p . The aim of the present section is 
to calculate the semiclassical form factor at small arguments (31), using the 
result (26) for the distribution of pencils of periodic orbits in the barrier 
billiard. Let us take a C°°, compactly supported test function and integrate 
the distribution K di ^(r) over r. If the density of periodic pencils depends 
linearly on / (as is the case for the barrier billiard or the rectangular billiard), 
the integration over families of periodic orbits yields i^ diag (r) = A@(r), where 
A is a constant and the Heaviside step function. In such a case, we define 
K 2 (0) = A. 

As an introduction, we first deal with the simpler case of a rectangular 
billiard. 



12 



4.2 Rectangular billiard 



In the case of the rectangular billiard, discussed in section 2.1, the periodic 
orbits have lengths nl pp , where l pp is given by (1) with (M,N) coprime, 
and n e N is the repetition number. The area of each pencil of primitive 
periodic orbits pp is A pp = AA. When the sides a and b of the rectangle 
are incommensurable, there is only one pencil of length l pp and therefore in 
Equation (31) S p = AA. Equation (31) becomes 

^w = ^£E^V. -***/»> (32) 

pp n Fl 

hence (using the fact that d = A/ An and turning the sum over pp = (M, N) 
with M and N coprime into an integral over / with density p pp (l)) 

Q A POO 1 

K di ^(r) = —X^ dl —p (1)5(1 - Ankdr/n) (33) 

The density p pp (l) of periodic orbits is given by Equation (3) and yields 
-fT diag (r) = 1, as expected for integrable systems. 



4.3 Barrier billiard 

In the case of the barrier billiard, the periodic orbits have a length of the 
form nwl p with l p given by (5): here the primitive length is wl p and n is the 
repetition number. Two pencils of periodic orbits p and p' have the same 
length provided there exist repetition numbers n and n' such that nwl p = 
n'w'lpi. When a and b are incommensurable, this implies p = p', i.e. two 
pencils can have same length only if they are in the same direction. For a 
given direction (M, N) with M and iV coprime (which will now be labeled 
by p), there are three cylinders of area Ai and length Wil p , 1 < i < 3, and 
therefore w, w' belong to the set {wi,w 2 ,Wi + w 2 }- Equation (31) becomes 

^W = i:EE^K-«» (34) 

p n F 

where l p is given by (5) and S PjU is the sum over the Ai corresponding to a 
Wi which divides n: 



Sp t n — ^ ^ AjS Wi \ n , (35) 



i=l 
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with S r \t = 1 if r divides t, otherwise. Each area Ai is equal to (2b/M)hi x 
(wilp) cos ip p (ip p is the angle between the orbit and the horizontal). This can 
be rewritten as 

4A 

Ai = — hiWi (36) 
Q 

(note that since £V hiWi = s\Wi + s 2 w 2 = 9, one has £V = 4^4, i.e. 
the total area of the translation surface, as expected). Therefore S P:Tl only 
depends of the five numbers / = (wi, w 2 , si, s 2 ) and /13, and can be rewritten: 

AA 

Sf,h 3 ,n = [(Si - /l3)lWl4i|n + («2 ~ h)w 2 5 W2 \ n + /l 3 (u>i + W 2 )5( t01 + tl , 2 )|„] . 

q (37) 
The sum (34) over all periodic orbits can be partitioned into sums running 
over primitive pencils of periodic orbits p(f, h) belonging to a family / with 
a height of the long cylinder in [h, h + dh[; (34) becomes 

f&F J p(f,h) n p 

Each of the sums corresponding to a family / can be replaced, as in (33), by 
an integral with density p pp ,f,h(l), and S PjTl by Sf^,n- 

n = g| ih p^£jy p^vi s(l _ i^i). (39) 

Replacing the density p pp j,h by its expression (26), the integration over / 
becomes straightforward and yields K diag (r) = K 2 (0)Q(t), where 

k*® = ^ E h £ / dh [{S1 " (40) 

+ (s 2 - h)w 2 S W2 \ n + h(Wi + W 2 )5(«, 1+ «, 2 )|n] 2 C/ )h 

(we have used the fact that d = A/ 'Att). Expanding the square, we can 
perform the summation over n, using the identity 

7T 2 gcd(Wi,W 2 ) 2 

^ n 2 6 w 2 w 2 

n=l 1 z 

The form factor can therefore be written, after simplifications using the fact 
that gcd(u>i, w 1 + w 2 ) = gcd(u? 2 , wi + w 2 ) = gcd(w 1 ,w 2 ), as 



T<m = \T,f dh <*> k - 2 f * + * + " gcd( r" T ] \ ) h 

(42) 



+ s ; + a i + 28lS2 gcd(wi ' W2)S 

W1W2 
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Replacing the weight c/^ by its expression (25), we can easily perform the 
integration over h, which consists of terms of the form 

p^) min(tl ,^ r . 

for z/ = 0, 1, 2. The form factor becomes 

K 2 (0) = ^Wi«J 2 (wi + w 2 ) [min(si,s 2 ) 3 

q (si."2) =1 r|(lUl,t«2) 
s 1 jii 1 +s 2 '"2=9 

f , , gcd(w 1 ,w 2 ) 2 \ . 2 
- Si + S2 + g 7 ; r mm(s 1 ,s 2 ) 

+ f 8 ; + s l + 2 8l82 gCd( ""-'" 2)2 )min( 8l ,, 2 )l. (44) 

This sum can be evaluated with some cumbersome arithmetic manipulations; 
the calculation is given in the Appendix, and the final result is unexpectedly 
simple: 

1 1 



K 2 (0) = - + - (45) 
z q 

There are several comments to make concerning this value. First, it is close 
to the result corresponding to semi-Poisson statistics K 2 (0) = 1/2 [12]. This 
result is not valid for q = 2, since in that case there is an additional sym- 
metry in the billiard, with respect to the barrier, and the spectrum has to 
be desymmetrized. The calculation in this case has been done in [29] for a 
height of the barrier equal to 6/2 (half the height of the rectangle), and yields 
K 2 (0) = 1/2. The calculation for q = 2 and a barrier with any height has 
been done in [18] using a different method, and also yields K 2 (0) = 1/2. 

The result (45) is similar to previously obtained results [12] for rational 
polygonal billiards having the Veech property. For instance for triangular 
billiards with angles (vr/2, ir/n, tt/2 — ir/n) the form factor at the origin was 
found to be between 1/3 and 3/5 [12]. Here the form factor lies between 1/2 
and 5/6, which again is close to the semi-Poisson result. 
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Appendix 

In this appendix, we want to evalute the quantity 

^) = 4r E E ^M^2,^i,^,?), (46) 

a 1 w 1 +s 2 m 2 = q 

where 

f(s 1 ,s 2 ,w 1 ,w 2 ,q) = w 1 w 2 (w 1 + w 2 ) [min(si,s 2 ) 3 

. gcd(wi,w 2 ) 2 \ . 2 
Si + S2 + 9 7 ; r mm(si, s 2 ) 

WiW 2 {Wi + W 2 J J 
( 2 , 2 , o gcd(Wl, W 2 f \ . ] 

+ s x + s 2 + 2sis 2 mm(si,s 2 ) .(47) 

V " WiW 2 J 

The function / is homogeneous, in the sense that it verifies 

f(si, s 2 , Xw!, \w 2 , Xq) = f(Xsi, Xs 2 , W!,w 2 , Xq) = A 3 /(si, s 2 , w ± ,w 2 , q). (48) 

In (46), the first sum goes over all integers Wi > 1 and Sj > 1, % — 1,2, 
verifying S\W\ + 52^2 = q an d gcd(s 1? s 2 ) = 1. The number iV 9 is given by 
(14). The first step is to exchange the sum over (s^iUj) and the sum over r 
in (46), and substitute Wi = ru{. using the homogeneity of /, we get 

KM=\Y.^y E /(*i»*2,«i,«2,;), (49) 

q r\q («1.S2) =1 

s 1 u 1 +s 2 u 2 =9/r 

To get rid of the co-primality condition on (si,S2) we use the exclusion- 
inclusion principle, which for any function gives 

oo oo 

<p(s,s')= 5>(*M*M«')- (50) 

(s,s')=l s,s'=l *=1 

This allows to rewrite the form factor as 
1 

^(0) = ^rE^W^E^W E /(^l,^2,«l,«2,^) (51) 



i- ^/x(r)r 2 ^//(*)* 3 E /(si,*2,«i,' 



r|g t|q si«i+s 2 «2='?/(»*t) 



Here the sum over t from 1 to oo has been replaced by a sum over t\q since 
for all the other values of t there is no value of (s 1 , s 2 , u±, u 2 ) fulfilling the 
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condition ts\Ui + ts 2 u 2 — 0.1 r - Again, the homogeneity of / (Equation (47)) 
has been used. Setting d = rt we get 

^(o) = ^ E ( E MOM*)** 2 * J E «i, «2, |). 

9 d|g \ t|d / siui+S2U2=q/d 



We need to evaluate 



where 



**(0) = irE f E^OM^d 2 * I {G q/i + H q/i ), 

q q d\q \ t\d ' J 



G n = ^2 u 1 u 2 (u 1 + u 2 )[rmn(s 1 ,s 2 ) 3 

siUi+S2U2=n 

- (si + s 2 ) min(si, s 2 ) 2 + (s? + s 2 + 2siS 2 ) min(si, s 2 )] 



(52) 



(53) 



and 



H n = E M 1 M 2(«1 + M 2 ) 



siui+S2«2=" 



-9 



2sis 2 



gcd(M!,M 2 ) 2 
UiU 2 {ui + M 2 ) 

gcd(Mi,M 2 ^ 



min(si, s 2 ) z 



MlM 2 



min(si, s 2 ) 



(54) 



The quantities G n and H n will be evaluated separately. This evaluation will 
require the use of a theorem proved in [23]: 
Theorem. Let / : Z 4 — > C such that 

/(a, 6, x, y) - /(x, y, a, 6) = /(-a, -6, x, y) - /(x, y, -a, -6) (55) 

for all integers a, b, x and y. Then for n G N, n > 1, 

a,b,x ,y > 1 

- /(a, a + 6, y - x, y) + f(b - a, b,x,x + y)- f(a + b,b,x,x~ y)\ 

d-i 

= E E [/(°» 5' x? rf) + /( ^' °' x) + /( ^' 5' x ' - x) 



>,\n x=l 



T) T) 71 71 

- /(:r, x - d, ^) - /(x, d, 0, -) - /(d, x, ^, 0) . (56) 
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a. Evaluation of G n 

We can immediately point out that the identity 

min(s 1 , s 2 ) 2 - («i + s 2 ) min(si, s 2 ) = -SiS 2 , (57) 

valid for any integers s± and s 2 , allows to simplify G n . We now need to 
evaluate the sum 

G n = ^2 u i u 2(ui + u 2 ) min(si, s 2 ) (s? + s 2 ~ s i s 2) • (58) 

s\u\+S2U2=n 

for any integer n. Writing min(a, b) = k(a + b — \a — b\), we have 



SlMl+S2M2=n 



O. = 5 E 

- ("Ui + M 2 )(S^ + S 2 - SiS 2 )|"Ui - M 2 



s?w 2 + s 2 ui - ^(s?ui + s 2 u 2 ) (59) 



2 ^ 

+ - ^2 uiu 2 {s\ui + s\u 2 ). 



3 

S1U1+S2«2=™ 



The first sum in (59) can be evaluated by applying Theorem (56) to the 
function 

f(a, b,x,y) = Uxy-^]\(a- b) (x - y) | {a 2 + b 2 - ab) (60) 



3 

and is equal to 

» (61) 



n 2 (n - 1) 



18 

d\n 

The second sum in (59) can be evaluated by applying Theorem (56) to the 
function f(a,b,x,y) = b 2 y^ — b 2 xy z (see [23]). It gives 

\ E° 3 ^ = ^E( 3rf3 +( 1 - 4 ^)- (62) 

ax+by=n d\n 



Finally we get 



n 2 



6 

d\n 



If we now evaluate the quantity 

y~] UiU 2 (ui + u 2 ) min(si, s 2 ) = \- ^ u x u 2 



2 

siui+S2U2=n s\ui+S2U2=n 



S\U 2 + S 2 U\ 



-^( s l«l + S 2 U 2 ) - {Ui + u 2 )\ui - U 2 \ 



2n / <n 

+ -T- M i M 2, (64) 



3 

SlUl+S2«2=« 
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the first sum is given by Theorem (56) applied to the function 

1 / xy \ 

f(a,b,x,y) = - [xy — \(a - b)(x - y)\ (65) 



3 V 2 

and the second one is given by Theorem (56) applied to the function 
f(a,b,x,y) = nxy/3; altogether, this gives 

^ uiu 2 {ui + u 2 ) min(si, s 2 ) = ^ ^ (d 3 - nd) . (66) 

s 1 u 1 +s 2 U2=n d\n 

Together with Equation (63) we get 

Gn = \ ^2 u i u 2{ui + u 2 )mm(s 1 ,s 2 ). (67) 

s\u\+S2U2=n 

b. Evaluation of H n 

We want to evaluate 

H n = UiU 2 {ui + u 2 ) min(si, s 2 ) 

s\u\+S2U2=n 

min(s 1 ,s 2 ) 2s!S 2 \ 2 
-n ■ r + gcd ( Ul , u 2 ) 68 

UiU 2 {Ui + U 2 ) UiU 2 J 

for any integer n. Summing over all the possible values r of the gcd of u\ 
and u 2 , and substituting Ui = rvi, we have 

H n =y y r v^vx + v 2 ) mm(s 1 , s 2 ) ■ H . 

i , V ?"^2(Vl + V2) ^1^2/ 

(69) 

Then, as before, the co-primality condition can be expressed by a sum over 
t (see Equation (50)). Restricting the sum over t as before, we get 

H n = ^2^2fi(t)r 3 t ^ V!V 2 (vi + v 2 )mm(si,s 2 ) 

r\n t\n siVi+S2V2=n/(rt) 

n min(s 1 ,s 2 ) | 2s!S 2 \ 



rt viv 2 (vi + v 2 ) v±v 2 



Setting d = rt we get 



i 2 , 

!|n \ t|d / si<;i+S2f2=n/d 



n min(si,s 2 ) + 2_£i£ 2 \ 



dviv 2 {vi + i; 2 ) fiw 2 / ' 
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Let us now evaluate, for any integer m, the quantity 

v ST ( , ^ ■ ( \( min(si,s 2 ) 2s l s 2 \ 
K m = > ^1^2(^1 + v 2 ) mm(s 1 , s 2 ) -m ■ r H 

f— ' V ^1^2(^1 + ^2) ui«2 / 

Sll>l+S2«2=m 

^ min(si, s 2 ) 2 + 2 ^ sis 2 (ui + v 2 ) min(si, s 2 ). (72) 



-m 

Slt)l+S2D2=Tl Sl«l+S2f2=m 



Let 



^ (2si« 2 - viv 2 ) {vi + v 2 ) min(si, s 2 ) (73) 

51^1+52^2=™ 

= ^2 (2^1^2(51 + s 2 ) min^i, w 2 ) - ^i^^i + v 2 ) min(si, s 2 )) 

Slfl+S2^ , 2=™ 

after exchanging (s 1; s 2 ) and (fi,f 2 ) in the first half of the right member. 
Writing min(a, b) = \{a + b — \a — b\) and applying Theorem (56) to the 
function ^ 

f(a,b,x,y) = -- \ab(a - b)(x - y)\ (74) 

one gets 

d-l 

L m = m^2^2 x(d - x). (75) 

d\m x=l 

Applying Theorem (56) to the function 

f(a,b,x,y) = ^ (76) 

one gets 

d-l 

min(s 1 , s 2 f = x (d - x) . (77) 

s 1 v 1 +s 2 V2=m d\m x=l 

This proves that 

K m= ^ V\V2{v\ + V2) min(si, s 2 ) (78) 

s\v\+S2V2=m 

and therefore 

( d A 

ff » = E Yl ^i^K +^2)min(si,s 2 ). (79) 

d\n \ t\d I s±vi+S2V2=n/d 
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c. Calculation of ^(0) 

The evaluation of (52) will require to introduce the functions 

/(„) = **!> and 9W = E# m 
n d z 

d\n 

For fi and / 2 two arithmetic functions, the Dirichlet convolution is defined 
by 

/i*/2(n) = £/i(d)/ 2 (^). (81) 
Replacing the expressions found for G q / d and H q / d in Equation (52) we get 

Uiu 2 {ui + u 2 ) min(si, s 2 ) 



+ ^ d' 3 g(d') uiu 2 (ui + U2) min(si, s 2 ) 

d'\q/d siU\+S2U2=q/dd' 



) ■ (82) 



Rewriting the constant N q given by (14), using the inclusion-exclusion prin- 
ciple and following the steps from Equations (49) to (52), we get 

N q = J^(/ * n)(d)d 2 ^ «iW2(wi + «2)nun(si,s 2 )- (83) 

d\q siUi+S2U2=q/d 

We see that the first term in (82) is equal to 1/2. If we set 5 = dd', the 
second term gives 

E^) 3 £(/ * v)(d)g( 6 j) u ' u ^ + M2 ) min ( Sl ' S2 )- ( 84 ) 

' <5|g \ d|<5 / si«i+s2M2=g/<5 

But 

£(/ * ^)^) = K/ * * = [/ * * (85) 

d|<J 

by associativity of Dirichlet convolution, and 

(/**</)(*) = (86) 

by Moebius inversion formula. Finally the term (84) simplifies to 1/q, which 
completes the proof. 
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